Low-noise conditional operation of singlet-triplet coupled quantum dot qubits 
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We theoretically study the influence of charge noise on a controlled phase gate, implemented using 
two proximal double quantum dots coupled electrostatically. Using the configuration interaction 
method, we present a full description of the conditional control scheme and quantitatively calculate 
the gate error arising from charge fluctuations. Our key finding is that the existence of noise-immune 
sweet spots depends on not only the energy detuning but also the device geometry. The conditions 
for sweet spots with minimal charge noise are predicted analytically and verified numerically. Going 
beyond the simple sweet-spot concept we demonstrate the existence of other optimal situations for 
fast and low- noise singlet- triplet two-qubit gates. 
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Coupled quantum dots are promising candidates for 
future implementations of quantum computatioi^i — . 
They have potentially excellent scalability due to the 
well-developed semiconductor nanoelectronics technol- 
ogy. Moreover, the confinement potential can be electri- 
cally tuned by nearby lithographic gates, enabling easy 
controllability of quantum dynamics^^^. However, such 
electrostatic controllability also makes the system vul- 
nerable to electrical fluctuations in the environmenlii^rii, 
leading to decoherence, thus hindering the requisite co- 
herent manipulation of quantum states. To minimize the 
influence of charge noise, one must search for optimal 
conditions or "sweet spots" in the parameter spacei^ — , 
where the first-order response of exchange energy to 
voltage fluctuations is zero, so that the decoherence of 
the exchange-coupled spin qubits to the environmental 
charge noise is minimal. 

Recently, much attention has been focused on the 
multi-qubit coherent operations in coupled quantum 
dot^i^— . This is partly because all the needed single- 
qubit manipulations have already been demonstrated in 



experiments 



8-10 



and partly because implementing ele- 



mentary multi-qubit gates would be an important mile- 
stone en route to a scalable quantum computer. In a re- 
cent experiment, van Weperen et al?^ demonstrated the 
fast conditional operation of a singlet-triplet qubit con- 
trolled by an adjacent two-electron double quantum dot. 
A two-qubit controlled-phase (C-phase) gate is realized 
through the capacitive interaction between the double 
quantum dots (DQDs). However, the problem of charge 
fluctuations in this scheme is severe. It is therefore of vi- 
tal importance to quantitatively understand charge noise 
and find multiqubit "sweet spots" for improving the per- 
formance and reliability of the two-qubit gate. 

The goal of this theoretical work is to find sweet spots 
for the two-qubit conditional operation of the so-called 
singlet-triplet qubits^^^^^^^, so that the impact of charge 
fluctuations is minimized. In particular, we explore for 
the first time the significance of the device geometry in 
reducing charge noise. In comparison to the previous 
investigation^—, a general arrangement of quantum 
dots with arbitrary geometry and relative configuration 
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FIG. 1. (Color online) (a) General configuration of the con- 
trol and target double quantum dots, (b) Loading electrons 
into the control DQD effectively changes the bias and tunnel 
coupling of the target DQD. (c)-(d) When the control DQD 
is tuned from (0,2) to (1,1) charge states, the variations of 
As and At are plotted as functions of Oc and Wr- 



angles [Fig.[Ha)] is studied in detail. Based on the multi- 
electron configuration interaction calculations of the cou- 
pled double-dot system, we address the following ques- 
tions: (i) How does the control mechanism depend on 
the device geometry? (ii) How is the C-phase gate af- 
fected by charge noise? (iii) Which configurations sup- 
port sweet spots? (iv) What are the optimal detuning 
energies? (v) Are there any other favorable cases? The 
answers to these questions could be of considerable help 
to the ongoing double-dot multiqubit experiments in var- 
ious laboratories worldwide. 

We first introduce a microscopic model and provide 
a clear understanding of the conditional control scheme. 
In general, the Hamiltonian of a quantum dot system is 
given by i:^ = Y.ki ^kiclci + Y.hjki Ghjkicic]ckCu where 
the sum includes all possible terms that conserve the to- 
tal particle number and the total spin^s^. In our case, 



quantum tunnel couplings between the control and tar- 
get DQDs are highly suppressed because of the large po- 
tential barrier between the two qubits with each double- 
dot system being the usual singlet-triplet qubit, so that 
they are only capacitively coupled^^. Hence, we have 
a simplified Hamiltonian H = He + Ht + i^int- Here 
He and Ht describe the isolated control and target 
DQDs^— , and Hint describes their mutual interactions. 
More specifically, Hi^t includes three parts: the classi- 
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the occupation- modulated hoppings Hjt = — (J^i^i + 

Jt2n2) Ea(4aC4a + H.C.) - ( J^3^3 + -^^^4) Y.^{c\^C2a + 

H.C.), and the hopping- hopping interactions Hjh = —Jh 
^i 3 i' 3' ^i^j^i'^j' •> where i (j) and i' (/) are chosen from 
{1,2} '({3,4}) and i j^ i' {j j^ f). The Hjh term is ig- 
nored hereafter since it plays the same role as Hjt on the 
target qubit, but is two orders of magnitude smaller. 

Therefore, as shown in Fig. [T]^b), the control DQD in- 
fluences the target DQD in two ways: (i) changing the 
energy detuning. 



Ae^i'i) = U23 + U24 - Uis - t/i4, 

^^(0,2) _ 2/724 - 2/7i4, 



and (ii) changing the tunnel coupling, 

Ai(l'l) = Jt3 + Jt4, At(0'2) = 2Jt4. 
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Here the superscripts (1, 1) and (0,2) denote the charge 
states of the control DQD, which are determined by its 
spin state (triplet Tq or singlet S) via Pauli blockade^i^. 
Taking Eqs. ([T]) and ([2]) into account, the Hamiltonian of 
the target DQD is rewritten in a matrix form 



(3) 



where the basis sets are {C'l^Cj, |0), [{Cl ^Cl \ -\- 

^2,t^l,i) /v 2] 10)7 ^2,t^2,i|0)' [(^l,t^2,i ~ ^2,t^l,i) 

/\/2] |0)}, U is the on-site Coulomb interaction, V is 
the inter-site Coulomb interaction, e^^ = £j^ -\- As", 
^eff = ^ + ^^ + ^^"' ^ = (1' 1) o^ (0^2), and the defi- 
nition of Je, Jpj Jtj and St are shown in footnote |27| . 
Diagonalizing the above Hamiltonian gives the exchange 
energy J" , which is the energy difference between the two 
lowest eigenvalues (one singlet and the other triplet). As 
a result, the coherent precession of the target qubit is 
controlled by the charge state of the control DQD via 

Using the configuration interaction method, all the 
coupHng parameters in Eqs. ([I])-© can be readily calcu- 
lated for a given confinement potential using the lowest- 
energy Fock-Darwin stdites^^^^^^. In this paper, we 
adopt the quadratic confinement potential V{x^ y) = min 
[W{xi,yi) -111,'-' ,W(x4,^4) -^4,0], where W{xi,yi) 
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FIG. 2. (Color online) (a)-(b) Two typical configurations of 
the control and target DQDs. (c)-(d) Differences between the 
exchange energies with respect to the bias of the target DQD 
St- (e)-(f) Relative errors of the controlled phase gate with 
respect to £t- (g)-(h) Close-ups of the low- noise spots. The 
panels on the left are for configuration (a), those on the right 
are for configuration (b). 



= mujQ [{x — Xi)'^ -\- {y — ^i)^]/2 represents the quantum 
dot centered at (xi^yi). The numerical calculations are 
carried out for a GaAs systen>iiii^ with hwo = 3.96 meV, 
a = 30 nm, and d = 100 nm [Fig. [I{a)]. 

In Fig. [Hc)-[I];d), the results of Ae^^'^) - Ae^^^^^ and 
At^^'^) - At(^'2) are plotted as functions of Oc and Ot- 
Their different patterns indicate that the ratio R = 
[Ae^i'i) - A£:(0'2)]/[At^i'i) - At^^'^)] varies with the de- 
vice geometry. Therefore, the geometric configuration of 
the system determines whether the energy detuning or 
the tunnel coupling plays a more important role. As an 
example, we consider two typical configurations shown in 
Fig.[21^a)-|2];b). For configuration (a) with Oc = Ot = 0, 
Fig. mc)-[Hd) gives i? ^ 8, which means that the energy 
detuning dominates the control process. However, for 
configuration (b) with Oc = ^ and Ot = 7r/2, one finds 
i? = 0, indicating that the tunnel coupling plays the de- 
cisive role. As will be shown later, the ratio R is crucial 
for reducing charge noise. 

Next, we quantitatively investigate the gate error aris- 
ing from charge noise. To facilitate the following discus- 
sion, we define a Bloch sphere for the target qubit, with 
l^*) being the north pole and \Tq) being the south pole^i^i. 
Following the experiments^, the coherent manipulation of 
the target qubit consists of the following three steps: (i) 



preparing an initial state and adiabatically loading it into 
the x-y plane of the Bloch sphere, (ii) rotating it about 
the z-axis through an angle 0^ = J^r/h during the pre- 
cession time r, and (iii) adiabatically unloading it out 
of the x-y plane and measuring the final state. We note 
that in step (ii), different control states result in differ- 
ent rotation angles 0^. This difference [0 = O'^^'^^ _ 5/(0)2)j 
finally gives rise to a controlled ^-phase gate. 

In reality, unavoidable environmental charge noise af- 
fects the confinement potential of quantum dots, and 
therefore perturbs the exchange energy J". For ex- 
change errors 6J^^'^^ and SJ^^''^\ becomes 0-\-S0^ where 
60= [(^J^^'^^— ^J^^'^^]r/fi. If we aim to perform a ^-phase 
gate, the gate time should be r = M/| J*^^'^^ — J*^^'^^|, and 
the relative error is given by 

r, = \S0/e\ = |,5j(i-i) - ,5J(°-2)|/| J(i'i) - J(0'2)|. (4) 

Although there are many sources of charge 
noise^^^^i^ (background charge noised— i^, 
gate noiseii, etc.), their effects on the target DQD are 
similar: raising or lowering the central barrier Vb (thus 
changing in tunnel coupling t) and detuning the energy 
difference St- Thus we have 
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(5) 



77 = defines a sweet spot for the target qubit. To eval- 
uate 77, we make the reasonable assumption that SVb 
and Sst are of the same magnitude and independent^^. 
For GaAs quantum dots the fluctuation is found to be 
SVb ^ Set ^ 0.07 - O.W/j.eY'^^^. We choose SVb = 
Sst = 0.2/ieV in our calculation, which corresponds to 
the worst case scenario with stronger fluctuations. All 
the other terms in Eq. (J5]) can be readily calculated us- 
ing the configuration interaction method^i^i^i^i^. 

So far we have shown how charge fluctuations affect 
the conditional operation. In the following, we will find 
optimal conditions, or "sweet spots", that suppress r]. 
For practical purposes, it would be advantageous to keep 
I j(i,i) _ j(0'2)| appreciable at the same time. This turns 
out to be possible as we show below. In Fig. [2lc)-|2ff), 
we plot calculated J^^'^^ — j(^'^) and r] as functions of 
St for the two typical configurations. For configuration 
(b), there exist two sweet spots near the crossovers of 
charge sectors. At the sweet spots, the relative error r] 
becomes zero, and |J^^'^^ — J^^'^^| reaches its maximum 
value. This means that if Et can be fixed precisely at the 
sweet spots, we get a robust and relatively fast C-phase 
gate. In panel (f), 77 seems to be sensitive to St in the 
vicinity of sweet spots. However, its close-up [panel (h)] 
shows that 77<5xlO~^as long as the error of St is no 
more than ±0.08 meV. Thus the low-noise gate can be 
performed with the state-of-the-art techniques. On the 
contrary, configuration (a) does not support any sweet 
spot in the (1, 1)t charge sector. Only an optimal point 
with a relatively small rj is found at the blue square in 
panel (e). However, configuration (a) also offers some 
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FIG. 3. (Color online) Low-energy spectra of the target DQD 
near the crossovers of the (2,0)t and (1, 1)t charge sectors 
[panel (a)], and the (1, 1)t and (0, 2)t charge sectors [panel 
(b)]. The superscripts denote the charge states of the control 
DQD. Insets: The corresponding J^^'^^ — J*^°'^) with respect 
to St- Two sweet spots are found at st = s^ ^.^. The dashed 
lines denote the boundaries between different charge sectors. 
A sweet spot in (1, 1)t charge sector is good for reducing 
dephasing error. 



advantages. First, panels (c)-(d) show that the value of 
|j(i,i) _ j(0'2)| -^^ (.g^gg (^g^j jg usually larger than the one 
in case (b), leading to a faster gate operation. Second, 
as can be seen from panel (g), working at the optimal 
point does not require a precise control of Et- In ad- 
dition, the minimal possible r] (r^min = 1.36 x 10~^) is 
rather small, allowing in principle the possibility of fault- 
tolerant quantum computation. To summarize, we have 
shown the existence of sweet spots in case (a). We have 
also gone beyond the sweet spots and found an optimal 
point of operation for configuration (b). 

Now we turn to the general configurations and provide 
analytic expressions for the sweet spots. The low-energy 
spectrum of the target DQD is shown in Fig. [3Ka)-[3l^b). 
We note that near the crossovers of different charge sec- 
tors, the singlet states can be described approximately 
by a two-level anti-crossing model 
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where A^ = {U — V ^ 6:^^)/2, the "+" superscript de- 
notes the (2,0)t and (1, 1)t crossover [see panel (a)], 
and the "-" superscript denotes the (1,1)t and (0, 2)t 
crossover [see panel (b)]. Here we have neglected the 

small Je and Jp terms for simplicity. Diagonalizing Hr^' 
gives the energy of the lowest singlet state E^^^ = A J 
— Y^2(tg^)^ + (^±)^- On the other hand, as shown by the 
horizontal lines in Fig. [3l the lowest triplet state has a 
constant energy. Since the exchange energy J" is the 
gap between the lowest singlet and the triplet state, we 
obtain [J(i'i) - J(0^2)]± _ ^(o,2),± _ ^(i,i),±. Accord- 
ing to Eq. dSJ, we can find the sweet spots by solving 

simultaneously. The sweet spots are given by 

4,min = MU -V)- Ae^i'i) + [t + Jt+ Ai(i'i)]i? (7) 
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FIG. 4. (Color online) (a) In the (1, 1)t charge sector, the 
minimal possible relative error ?7min as a function of Oc and 
Ot- (b) The optimal detuning of the target DQD with respect 
to Oc and Ot- 



where R= [Ae^^^^^ -Ae^^^^^]/[At^^^^^ -At^^^^^]. However, 
there are some caveats here. First, Eq. (J6j) is a good 
approximation only when the sweet spots are close to the 
crossovers ±{U — V) — Ae^ (dashed lines in Fig. [3|), i.e., 
the absolute value of R is small [see Eq. ([7])]. Second, we 
want to find the sweet spots in the (1, 1)t charge sector, 
i.e., V -U - min[A£(i'i), A£(0'2)] < £T,min <U-V- 
max[A£^^'^\ Ae^^'^^], because a high double-occupation 
probability would result in a large dephasing errori2i. 
In general, as can be seen from Eq. ([7]) and the insets of 
Fig. [31 only one sweet spot meets this requirement, given 
by the condition that R < (i? > 0), ^Jmin (^Tmin) ^^ 
inside the (1, 1)t regime. 

The above results are further verified by exact numer- 
ical calculations. As shown in Fig. HJ^a), the minimal 
possible relative error in the (1, 1)t charge sector is plot- 
ted as a function of Oc and Ot- Comparing the white 
regions with Fig. Hl^c), one finds that the system is im- 
mune to charge noise when Ae^^'^^ — Ae^^'^^ ^ (i.e.. 



R^ 0) and Et = ^T,min- In contrast, from the dark blue 
regions and Fig. HJc), we see that the system is always 
sensitive to charge noise when At^^'^"^ — At^^'^^ ^ (i.e., 
R ^ oo). Moreover, there are some light blue regions in 
between, where the system has an optimal point as the 
one shown in Fig. [2fa). In Fig. Hl^b), £:T,min of the sweet 
spots is also in excellent agreement with the analytical 
approximation in Eq. ([7]). Therefore, the ratio R and the 
device geometry together determine the sweet spots and 
their applicability in coherent qubit manipulations. 

In summary, we have studied the coherent multiqubit 
operations in the coupled DQD system with a general 
geometry emphasizing how to reduce charge noise. We 
have developed a microscopic model to fully understand 
the mechanism of the conditional operation. We have 
shown that the exchange energy of the target qubit is 
affected by the control qubit through two channels: en- 
ergy detuning and tunnel coupling. In particular, the 
device geometry determines which one plays a domi- 
nant role. Using the configuration interaction method, 
we have evaluated the relative error of the C-phase gate 
arising from charge noise. We have demonstrated the 
existence of noise-immune sweet spots in some optimal 
configurations, where the tunnel coupling serves as the 
main control channel. On the contrary, one always sees 
large charge fluctuations in the (1,1)t charge sector if 
the tunnel coupling makes no contribution. We have fur- 
ther developed a two-level anti-crossing model which an- 
alytically describes the conditions for sweet spots. In 
addition, we have found some optimal points in the pa- 
rameter space, where the gate is fast and insensitive to 
charge noise. Our work should guide future experimen- 
tal efforts to carry out coherent multiqubit operations in 
double-quantum-dot structures. 
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